Reinfection rate and disease severity of the BA.5 Omicron SARS-CoV-2 lineage compared to previously circulating variants of concern in the Canary Islands (Spain).

ABSTRACTThe emergence of the Omicron SARS-CoV-2 variant of concern has changed the COVID-19 scenario as this variant is characterized by high transmissibility and immune evasion ability. To evaluate the impact of this variant on the Canary Islands (Spain) population, we determined the reinfection rates and disease severity associated with the Omicron sublineages and the previously circulating variants of concern. We performed a retrospective observational study on 21,745 SARS-CoV-2 viral genomes collected from December 2020 to July 2022 in the Canary Islands (Spain). We compared the reinfection rates between lineages using pairwise proportion and Fisher's exact tests. To assess disease severity, we studied the association of Alpha, Delta, BA.1, BA.2, BA.5, and other risk factors on 28-day hospital mortality using logistic regression and Cox proportional hazard models. We observed 127 bona fide reinfection cases throughout the study period. We found that BA.5 had the highest reinfection rate compared to other lineages (vs. Delta p = 2.89 × 10-25; vs. BA.1 p = 5.17 × 10-11; vs. BA.2 p = 0.002). Among the 1,094 hospitalized patients, multivariate logistic regression showed that Alpha (Odds Ratio [OR] = 0.45, 95% Confidence Interval [CI] = 0.23-0.87, p = 0.02), BA.2 (OR = 0.38, 95% CI = 0.22-0.63, p = 1.91 × 10-4), and BA.5 (OR = 0.30, 95% CI = 0.16-0.55, p = 1.05 × 10-4) had lower 28-day hospital mortality compared to Delta. These results were confirmed by using Cox proportional hazard models. Omicron lineages, and in particular BA.5, were associated with higher reinfection rates and lower disease severity (28-day hospital mortality) than previously circulating variants of concern.

Bioinformatic approaches to draft the viral genome sequence of Canary Islands cases related to the multicountry mpox virus 2022-outbreak.

On July 23, 2022, monkeypox disease (mpox) was declared a Public Emergency of International Concern (PHEIC) by the World Health Organization (WHO) due to a multicountry outbreak. In Europe, several cases of mpox virus (MPXV) infection related to this outbreak were detected in the Canary Islands (Spain). Here we describe the combination of viral DNA sequencing and bioinformatic approaches, including methods for de novo genome assembly and short- and long-read technologies, used to reconstruct the first MPXV genome isolated in the Canary Islands on the 31st of May 2022 from a male adult patient with mild symptoms. The same sequencing and bioinformatic approaches were then validated with three other positive cases of MPXV infection from the same mpox outbreak. We obtained the best results using a reference-based approach with short reads, evidencing 46-79 nucleotide variants against viral sequences from the 2018-2019 mpox outbreak and placing the viral sequences in the new B.1 sublineage of clade IIb of the MPXV classification. This study of MPXV demonstrates the potential of metagenomics sequencing for rapid and precise pathogen identification.

NanoRTax, a real-time pipeline for taxonomic and diversity analysis of nanopore 16S rRNA amplicon sequencing data.

Background: The study of microbial communities and their applications have been leveraged by advances in sequencing techniques and bioinformatics tools. The Oxford Nanopore Technologies long-read sequencing by nanopores provides a portable and cost-efficient platform for sequencing assays. While this opens the possibility of sequencing applications outside specialized environments and real-time analysis of data, complementing the existing efficient library preparation protocols with streamlined bioinformatic workflows is required. Results: Here we present NanoRTax, a Nextflow pipeline for nanopore 16S rRNA gene amplicon data that features state-of-the-art taxonomic classification tools and real-time capability. The pipeline is paired with a web-based visual interface to enable user-friendly inspections of the experiment in progress. NanoRTax workflow and a simulated real-time analysis were used to validate the prediction of adult Intensive Care Unit patient mortality based on full-length 16S rRNA sequencing data from respiratory microbiome samples. Conclusions: This constitutes a proof-of-concept simulation study of how real-time bioinformatic workflows could be used to shorten the turnaround times in critical care settings and provides an instrument for future research on early-response strategies for sepsis.

Evaluation of a whole-exome sequencing pipeline and benchmarking of causal germline variant prioritizers.

Most causal variants of Mendelian diseases are exonic. Whole-exome sequencing (WES) has become the diagnostic gold standard, but causative variant prioritization constitutes a bottleneck. Here we assessed an in-house sample-to-sequence pipeline and benchmarked free prioritization tools for germline causal variants from WES data. WES of 61 unselected patients with a known genetic disease cause was obtained. Variant prioritizations were performed by diverse tools and recorded to obtain a diagnostic yield when the causal variant was present in the first, fifth, and 10th top rankings. A fraction of causal variants was not captured by WES (8.2%) or did not pass the quality control criteria (13.1%). Most of the applications inspected were unavailable or had technical limitations, leaving nine tools for complete evaluation. Exomiser performed best in the top first rankings, while LIRICAL led in the top fifth rankings. Based on the more conservative top 10th rankings, Xrare had the highest diagnostic yield, followed by a three-way tie among Exomiser, LIRICAL, and PhenIX, then followed by AMELIE, TAPES, Phen-Gen,  AIVar, and VarNote-PAT. Xrare, Exomiser, LIRICAL, and PhenIX are the most efficient options for variant prioritization in real patient WES data.

Tracing the trajectories of SARS-CoV-2 variants of concern between December 2020 and September 2021 in the Canary Islands (Spain).

Several variants of concern (VOCs) explain most of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic waves in Europe. We aimed to dissect the spread of the SARS-CoV-2 VOCs in the Canary Islands (Spain) between December 2020 and September 2021 at a micro-geographical level. We sequenced the viral genome of 8,224 respiratory samples collected in the archipelago. We observed that Alpha (B.1.1.7) and Delta (B.1.617.2 and sublineages) were ubiquitously present in the islands, while Beta (B.1.351) and Gamma (P.1/P.1.1) had a heterogeneous distribution and were responsible for fewer and more controlled outbreaks. This work represents the largest effort for viral genomic surveillance in the Canary Islands so far, helping the public health bodies in decision-making throughout the pandemic.

Admixture Mapping of Sepsis in European Individuals With African Ancestries.

Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction. Although the ancestral genetic background is a relevant factor for sepsis susceptibility, there is a lack of studies using the genetic singularities of a recently admixed population to identify loci involved in sepsis susceptibility. Here we aimed to discover new sepsis loci by completing the first admixture mapping study of sepsis in Canary Islanders, leveraging their distinctive genetic makeup as a mixture of Europeans and African ancestries. We used a case-control approach and inferred local ancestry blocks from genome-wide data from 113,414 polymorphisms genotyped in 343 patients with sepsis and 410 unrelated controls, all ascertained for grandparental origin in the Canary Islands (Spain). Deviations in local ancestries between cases and controls were tested using logistic regressions, followed by fine-mapping analyses based on imputed genotypes, in silico functional assessments, and gene expression analysis centered on the region of interest. The admixture mapping analysis detected that local European ancestry in a locus spanning 1.2 megabases of chromosome 8p23.1 was associated with sepsis (lowest p = 1.37 × 10-4; Odds Ratio [OR] = 0.51; 95%CI = 0.40-0.66). Fine-mapping studies prioritized the variant rs13249564 within intron 1 of MFHAS1 gene associated with sepsis (p = 9.94 × 10-4; OR = 0.65; 95%CI = 0.50-0.84). Functional and gene expression analyses focused on 8p23.1 allowed us to identify alternative genes with possible biological plausibility such as defensins, which are well-known effector molecules of innate immunity. By completing the first admixture mapping study of sepsis, our results revealed a new genetic locus (8p23.1) harboring a number of genes with plausible implications in sepsis susceptibility.

Targeted analysis of genomic regions enriched in African ancestry reveals novel classical HLA alleles associated with asthma in Southwestern Europeans.

Despite asthma has a considerable genetic component, an important proportion of genetic risks remain unknown, especially for non-European populations. Canary Islanders have the largest African genetic ancestry observed among Southwestern Europeans and the highest asthma prevalence in Spain. Here we examined broad chromosomal regions previously associated with an excess of African genetic ancestry in Canary Islanders, with the aim of identifying novel risk variants associated with asthma susceptibility. In a two-stage cases-control study, we revealed a variant within HLA-DQB1 significantly associated with asthma risk (rs1049213, meta-analysis p = 1.30 × 10-7, OR [95% CI] = 1.74 [1.41-2.13]) previously associated with asthma and broad allergic phenotype. Subsequent fine-mapping analyses of classical HLA alleles revealed a novel allele significantly associated with asthma protection (HLA-DQA1*01:02, meta-analysis p = 3.98 × 10-4, OR [95% CI] = 0.64 [0.50-0.82]) that had been linked to infectious and autoimmune diseases, and peanut allergy. HLA haplotype analyses revealed a novel haplotype DQA1*01:02-DQB1*06:04 conferring asthma protection (meta-analysis p = 4.71 × 10-4, OR [95% CI] = 0.47 [0.29- 0.73]).

Whole-Blood Mitochondrial DNA Copies Are Associated With the Prognosis of Acute Respiratory Distress Syndrome After Sepsis.

Acute respiratory distress syndrome (ARDS) is an inflammatory process of the lungs that develops primarily in response to pulmonary or systemic sepsis, resulting in a disproportionate death toll in intensive care units (ICUs). Given its role as a critical activator of the inflammatory and innate immune responses, previous studies have reported that an increase of circulating cell-free mitochondrial DNA (mtDNA) is a biomarker for fatal outcome in the ICU. Here we analyzed the association of whole-blood mtDNA (wb-mtDNA) copies with 28-day survival from sepsis and sepsis-associated ARDS. We analyzed mtDNA data from 687 peripheral whole-blood samples within 24 h of sepsis diagnosis from unrelated Spanish patients with sepsis (264 with ARDS) included in the GEN-SEP study. The wb-mtDNA copies were obtained from the array intensities of selected probes, with 100% identity with mtDNA and with the largest number of mismatches with the nuclear sequences, and normalized across the individual-probe intensities. We used Cox regression models for testing the association with 28-day survival. We observed that wb-mtDNA copies were significantly associated with 28-day survival in ARDS patients (hazard ratio = 3.65, 95% confidence interval = 1.39-9.59, p = 0.009) but not in non-ARDS patients. Our findings support that wb-mtDNA copies at sepsis diagnosis could be considered an early prognostic biomarker in sepsis-associated ARDS patients. Future studies will be needed to evaluate the mechanistic links of this observation with the pathogenesis of ARDS.

Admixture mapping analysis reveals differential genetic ancestry associated with Chagas disease susceptibility in the Colombian population.

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP), and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the major histocompatibility complex region [Odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.66-0.83, lowest P-value = 4.53 × 10-8]. The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95% CI = 0.90-0.97, P-value = 3.54 × 10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the human leukocyte antigen genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease.

Nanopore sequencing and its application to the study of microbial communities.

Since its introduction, nanopore sequencing has enhanced our ability to study complex microbial samples through the possibility to sequence long reads in real time using inexpensive and portable technologies. The use of long reads has allowed to address several previously unsolved issues in the field, such as the resolution of complex genomic structures, and facilitated the access to metagenome assembled genomes (MAGs). Furthermore, the low cost and portability of platforms together with the development of rapid protocols and analysis pipelines have featured nanopore technology as an attractive and ever-growing tool for real-time in-field sequencing for environmental microbial analysis. This review provides an up-to-date summary of the experimental protocols and bioinformatic tools for the study of microbial communities using nanopore sequencing, highlighting the most important and recent research in the field with a major focus on infectious diseases. An overview of the main approaches including targeted and shotgun approaches, metatranscriptomics, epigenomics, and epitranscriptomics is provided, together with an outlook to the major challenges and perspectives over the use of this technology for microbial studies.

NanoCLUST: a species-level analysis of 16S rRNA nanopore sequencing data.

SUMMARY: NanoCLUST is an analysis pipeline for the classification of amplicon-based full-length 16S rRNA nanopore reads. It is characterized by an unsupervised read clustering step, based on Uniform Manifold Approximation and Projection (UMAP), followed by the construction of a polished read and subsequent Blast classification. Here, we demonstrate that NanoCLUST performs better than other state-of-the-art software in the characterization of two commercial mock communities, enabling accurate bacterial identification and abundance profile estimation at species-level resolution. AVAILABILITY AND IMPLEMENTATION: Source code, test data and documentation of NanoCLUST are freely available at https://github.com/genomicsITER/NanoCLUST under MIT License. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Admixture mapping of asthma in southwestern Europeans with North African ancestry influences.

The prevalence of asthma symptoms in Canary Islanders, a southwestern European population from Spain, is almost three times higher than the country average. Because the genetic risks identified so far explain <5% of asthma heritability, here we aimed to discover new asthma loci by completing the first admixture mapping study in Canary Islanders leveraging their distinctive genetic makeup, where significant northwest African influences coexist in the European genetic diversity landscape. A 2-stage study was conducted in 1,491 unrelated individuals self-declaring having a Canary Islands origin for the 4 grandparents. Local ancestry estimates were obtained for the shared positions with reference data from putative ancestral populations from Europe, North Africa, and sub-Saharan Africa. Case-control deviations in local ancestry were tested for each ancestry separately using logistic regressions adjusted for principal components, followed by fine-mapping analyses based on imputed genotypes and analyses of the likely deleterious exonic variants. The admixture mapping analysis of asthma detected that local North African ancestry in a locus spanning 365 kb of chromosome 16q23.3 was associated with asthma risk at study-wide significance [lowest P = 1.12 × 10-4; odds ratio (OR) = 2.05; 95% confidence interval (CI) = 1.41-3.00]. Fine-mapping studies identified a variant associated with asthma, and results were replicated in independent samples (rs3852738, OR = 1.34; 95% CI = 1.13-1.59, P = 7.58 × 10-4). Whole exome sequencing data from a subset of individuals revealed an enrichment of likely deleterious variants among asthma cases in 16q23.3, particularly in the phospholipase Cγ2 (PLCG2) gene (P = 3.67 × 10-4). By completing the first mapping study of asthma in admixed populations from Europe, our results revealed a new plausible asthma locus.

Resource availability drives bacterial succession during leaf-litter decomposition in a bromeliad ecosystem.

Despite the growing number of investigations on microbial succession during the last decade, most of our knowledge on primary succession of bacteria in natural environments comes from conceptual models and/or studies of chronosequences. Successional patterns of litter-degrading bacteria remain poorly documented, especially in undisturbed environments. Here we conducted an experiment with tank bromeliads as natural freshwater microcosms to assess major trends in bacterial succession on two leaf-litter species incubated with or without animal exclusion. We used amplicon sequencing and a co-occurrence network to assess changes in bacterial community structure according to treatments. Alpha-diversity and community complexity displayed the same trends regardless of the treatments, highlighting that primary succession of detrital-bacteria is subject to resource limitation and biological interactions, much like macro-organisms. Shifts in bacterial assemblages along the succession were characterized by an increase in uncharacterized taxa and potential N-fixing bacteria, the latter being involved in positive co-occurrence between taxa. These findings support the hypothesis of interdependence between taxa as a significant niche-based process shaping bacterial communities during the advanced stage of succession.

Allelopathic Effects of Native Versus Invasive Plants on One Major Invader.

Allelopathy is defined as the effects (stimulatory and inhibitory) of a plant on the development of neighboring plants through the release of secondary compounds. Autoallelophaty is the beneficial or harmful effect of a plant species on itself. The allelopathic potential belonging to a native species could induce a biotic resistance against invasive plants, whereas allelochemicals released by exotic species could favor the establishment of invasive species (invasional meltdown). The aim of our study was to examine the potential allelopathic effect of four plant species on the target species Ludwigia hexapetala using two experiments. In the first experiment, we tested the allelopathic effect of root and leaf leachates of the two congeneric exotic species Ludwigia hexapetala and Ludwigia peploides on L. hexapetala, while in the second experiment, we studied the allelopathic effect of root and leaf leachates of a sympatric exotic species Myriophyllum aquaticum and of one native species Mentha aquatica on L. hexapetala. We measured the stem length to calculate the relative growth rate and four physiological traits (nitrogen balance index and flavonol, chorophyll, anthocyanin indices) of the target plants on a weekly basis. At the end of the experiment, we determined the aboveground and belowground biomass. We also counted the number of lateral branches and measured their lengths. We found that the root leachates of L. peploides and of Myriophyllum aquaticum had stimulated the synthesis of flavonols of L. hexapetala. Leaf leachate of L. hexapetala also stimulated its own flavonol synthesis. Also, the root leachate of L. peploides had stimulated the total biomass and length of lateral branches of L. hexapetala, whereas the production of lateral branches had been stimulated by root leachates of both Ludwigia species and by leaf leachate of Myriophyllum aquaticum. The autoallelopathy of L. hexapetala could explain its invasiveness. Both leachates produced by Mentha aquatica had no effect on the physiological and morphological traits of the invasive L. hexapetala and indicated no biotic resistance in the recipient community. The two invasive plant species Myriophyllum aquaticum and L. peploides could favor the establishment of L. hexapetala. These results suggested an "invasional meltdown."

NanoDJ: a Dockerized Jupyter notebook for interactive Oxford Nanopore MinION sequence manipulation and genome assembly.

BACKGROUND: The Oxford Nanopore Technologies (ONT) MinION portable sequencer makes it possible to use cutting-edge genomic technologies in the field and the academic classroom. RESULTS: We present NanoDJ, a Jupyter notebook integration of tools for simplified manipulation and assembly of DNA sequences produced by ONT devices. It integrates basecalling, read trimming and quality control, simulation and plotting routines with a variety of widely used aligners and assemblers, including procedures for hybrid assembly. CONCLUSIONS: With the use of Jupyter-facilitated access to self-explanatory contents of applications and the interactive visualization of results, as well as by its distribution into a Docker software container, NanoDJ is aimed to simplify and make more reproducible ONT DNA sequence analysis. The NanoDJ package code, documentation and installation instructions are freely available at https://github.com/genomicsITER/NanoDJ .

Simulated drought regimes reveal community resilience and hydrological thresholds for altered decomposition.

Future climate scenarios forecast a 10-50% decline in rainfall in Eastern Amazonia. Altered precipitation patterns may change important ecosystem functions like decomposition through either changes in physical and chemical processes or shifts in the activity and/or composition of species. We experimentally manipulated hydroperiods (length of wet:dry cycles) in a tank bromeliad ecosystem to examine impacts on leaf litter decomposition. Gross loss of litter mass over 112 days was greatest in continuously submersed litter, lowest in continuously dry litter, and intermediate over a range of hydroperiods ranging from eight cycles of 7 wet:7 dry days to one cycle of 56 wet:56 dry days. The resilience of litter mass loss to hydroperiod length is due to a shift from biologically assisted decomposition (mostly microbial) at short wet:dry hydroperiods to physicochemical release of dissolved organic matter at longer wet:dry hydroperiods. Biologically assisted decomposition was maximized at wet:dry hydroperiods falling within the range of ambient conditions (12-22 consecutive dry days) but then declined under prolonged wet:dry hydroperiods (28 and 56 dry days. Fungal:bacterial ratios showed a similar pattern as biologically assisted decomposition to hydroperiod length. Our results suggest that microbial communities confer functional resilience to altered hydroperiod in tank bromeliad ecosystems. We predict a substantial decrease in biological activity relevant to decomposition under climate scenarios that increase consecutive dry days by 1.6- to 3.2-fold in our study area, whereas decreased frequency of dry periods will tend to increase the physicochemical component of decomposition.

Allelopathic Effect of the Invasive Ludwigia hexapetala on Growth of Three Macrophyte Species.

The release of allelochemicals by plants can affect the performance of other organisms positively or negatively. We tested the effects of aqueous extracts and leachates derived from the leaves and roots of the invasive water primrose (Ludwigia hexapetala) on one submerged native species - Ceratophyllum demersum, and two exotic species - the submerged Egeria densa and the emergent growth form of Myriophyllum aquaticum. The effect of the aqueous extracts and leachates of L. hexapetala on photosynthetic yield, growth (i.e., relative growth rate, leaf area), root length, and length of the lateral shoots of each species were analyzed in spring and in autumn. In autumn, an allelopathic effect was established on the traits of the three macrophytes species. The root extracts stimulated leaf area and the photosynthetic yield of C. demersum and of E. densa, whereas leaf treatments (leachates and extracts) and root leachate reduced the leaf area of M. aquaticum. The autumnal root leachate of L. hexapetala decreased the relative growth rate of C. demersum, whereas it had no effect on the two others plants. The root extract increased the length of lateral branches of M. aquaticum in autumn, suggesting a positive effect of L. hexapetala on the lateral growth of M. aquaticum. Three main allelochemicals were identified in leaves: quercitrin, prunin, myricitrin. The concentrations of these allelochemicals were greater in the leaf extract taken from L. hexapetala in autumn than in spring, and those found in the leaf leachates for both seasons. This assessment of autumnal allelopathy could help to explain the patterns of plant community succession in invaded areas.

The predatory behavior of the Neotropical social wasp Polybia rejecta.

We experimentally studied the predatory behavior of Polybia rejecta (Vespidae, Polistinae, Epiponini) towards 2-88 mm-long insects attracted to a UV light trap. Foragers, which began to hunt at 6:30, selected 4-14 mm-long prey insects. Prey detection by sight by hovering wasps was confirmed using decoys. After the wasps landed and walked along a sinuous path, prey were detected by contact or from a distance (1-3cm). This was followed by seizure, stinging (contrarily to most other known cases), prey manipulation and retrieval. Prey that flew off might be caught in flight. The prey load, representing 30.7% of a forager's weight, was optimized by capturing up to six small prey or two medium-sized prey successively (both of which might be consumed in situ). The foragers cut off the wings of larger prey or cut them into two pieces and returned to gather the second piece. The handling time increased exponentially with the weight of the prey. Partial loading (i.e., retrieving a load much inferior to the maximum possible) was likely related to social facilitation, a form of nest-based recruitment that was demonstrated through the experimental elimination of local enhancement by removing foragers (both mechanisms favor the exploitation of favorable patches).